home *** CD-ROM | disk | FTP | other *** search
/ Mac-Source 1994 July / Mac-Source_July_1994.iso / C and C++ / Libraries / stringsearch / bmsource / lc.fwd.inc.c-A < prev    next >
Encoding:
Text File  |  1991-05-06  |  3.0 KB  |  147 lines  |  [TEXT/MPS ]
  1. /*
  2.     search routine generated by gen.
  3.     skip=lc, match=fwd, shift=inc
  4. */
  5. #include    "freq.h"
  6. #include    "sd.h"
  7. /*
  8.  * The authors of this software are Andrew Hume and Daniel Sunday.
  9.  * 
  10.  * Copyright (c) 1991 by AT&T and Daniel Sunday.
  11.  * 
  12.  * Permission to use, copy, modify, and distribute this software for any
  13.  * purpose without fee is hereby granted, provided that this entire notice
  14.  * is included in all copies of any software which is or includes a copy
  15.  * or modification of this software and in all copies of the supporting
  16.  * documentation for such software.
  17.  * 
  18.  * THIS SOFTWARE IS BEING PROVIDED "AS IS", WITHOUT ANY EXPRESS OR IMPLIED
  19.  * WARRANTY.  IN PARTICULAR, NEITHER THE AUTHORS NOR AT&T MAKE ANY
  20.  * REPRESENTATION OR WARRANTY OF ANY KIND CONCERNING THE MERCHANTABILITY
  21.  * OF THIS SOFTWARE OR ITS FITNESS FOR ANY PARTICULAR PURPOSE.
  22.  */
  23.  
  24. #ifndef    CHARTYPE
  25. #define    CHARTYPE    unsigned char
  26. #endif
  27. #define    MAXPAT    256
  28.  
  29. #include    "stats.h"
  30.  
  31. #ifndef    TABTYPE
  32. #define    TABTYPE    long
  33. #endif
  34. typedef TABTYPE Tab;
  35.  
  36. static struct
  37. {
  38.     int patlen;
  39.     CHARTYPE pat[MAXPAT];
  40.     Tab delta[256];
  41.     int loopoffset;    /* patlen-1 - skip loop char */
  42. } pat;
  43.  
  44. prep(base, m)
  45.     CHARTYPE *base;
  46.     register m;
  47. {
  48.     CHARTYPE *skipc;
  49.     register CHARTYPE *pe, *pb;
  50.     register int j, r;
  51.     register Tab *d;
  52.     double tmax, fr;
  53.     extern double tcmp;
  54.     extern option;
  55.  
  56.     pat.patlen = m;
  57.     if(m > MAXPAT)
  58.         abort();
  59.     memcpy(pat.pat, base, m);
  60.     skipc = 0;
  61.     stats.len = m;
  62.     if(sd[m] == 0){
  63.         for(j = m; sd[j] == 0; j--)
  64.             sd[j] = 12;
  65.     }
  66.     d = pat.delta;
  67.     tmax = 2+tcmp;
  68.     for(pb = pat.pat, pe = pb+pat.patlen-1; pb <= pe; pb++){
  69.         fr = (1+tcmp*freq[*pb])/sd[pb-pat.pat];
  70. #ifdef    STATS
  71. /*printf("i=%d: f=%.4f freq[%c]=%.5f sd=%.1f\n", (pb-pat.pat), fr, *pb, freq[*pb], sd[pb-pat.pat]);/**/
  72. #endif
  73.         if(tmax > fr){
  74.             tmax = fr;
  75.             r = pb-pat.pat;
  76.         }
  77.     }
  78.     pat.loopoffset = m-1-r;
  79. if(option >= 0) pat.loopoffset=option;
  80.     skipc = &pat.pat[m-1-pat.loopoffset];
  81.     for(j = 0; j < 256; j++)
  82.         d[j] = m-pat.loopoffset;
  83.     for(pb = pat.pat, pe = pb+m-1-pat.loopoffset; pb <= pe; pb++){
  84.         d[*pb] = pe-pb;
  85.     }
  86. }
  87.  
  88. exec(base, n)
  89.     CHARTYPE *base;
  90. {
  91.     int nmatch = 0;
  92.     register CHARTYPE *e, *s;
  93.     register Tab *d0 = pat.delta;
  94.     register k, s_offset;
  95.     register CHARTYPE *p, *q;
  96.     register CHARTYPE *ep;
  97.  
  98.     k = pat.patlen-1-pat.loopoffset;    /* k is char we loop on */
  99.     s = base+k;
  100.     e = base+n;
  101.     memset(e, pat.pat[k], pat.patlen);
  102.     s_offset = -k;
  103.     ep = pat.pat + pat.patlen;
  104.     while(s < e){
  105. #ifdef    STATS
  106.         k = d0[*s];
  107.         stats.jump++;
  108.         while(k){
  109.             stats.jump++; stats.step[k]++;
  110.             k = d0[*(s += k)];
  111.             stats.jump++; stats.step[k]++;
  112.             k = d0[*(s += k)];
  113.             stats.jump++; stats.step[k]++;
  114.             k = d0[*(s += k)];
  115.         }
  116.         if(s >= e)
  117.             return(nmatch);
  118. #else
  119.         k = d0[*s];
  120.         while(k){
  121.             k = d0[*(s += k)];
  122.             k = d0[*(s += k)];
  123.             k = d0[*(s += k)];
  124.         }
  125.         if(s >= e)
  126.             return(nmatch);
  127. #endif
  128. #ifdef    STATS
  129.         stats.slow++;
  130. #endif
  131.         for(p = pat.pat, q = s+s_offset; p < ep; ){
  132. #ifdef    STATS
  133.             stats.cmp++;
  134. #endif
  135.             if(*q++ != *p++)
  136.                 goto mismatch;
  137.         }
  138.         nmatch++;
  139.     mismatch:
  140.         s++;
  141. #ifdef    STATS
  142.         stats.step[1]++;
  143. #endif
  144.     }
  145.     return(nmatch);
  146. }
  147.